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The ability of the Smaller Alignment Index (SALI) to distinguish chaotic from ordered 
motion, has been demonstrated recently in several publications^^ Basically it is observed 
that in chaotic regions the SALI goes to zero very rapidly, while it fluctuates around a 
nonzero value in ordered regions. In this paper, we make a first step forward explaining 
these results by studying in detail the evolution of small deviations from regular orbits lying 
on the invariant tori of an integrable 2D Hamiltonian system. We show that, in general, 
any two initial deviation vectors will eventually fall on the "tangent space" of the torus, 
pointing in different directions due to the different dynamics of the 2 integrals of motion, 
which means that the SALI (or the smaller angle between these vectors) will oscillate away 
^ from zero for all time. 

in 

m 

■ §1. Introduction 

| The evaluation of the Smaller Alignment Index (SALI) is an efficient and 

simple method to determine the ordered or chaotic nature of orbits in dynamical 
systems. The SALI was proposed in Ref. ^) and it has been successfully applied to 
distinguish between ordered and chaotic motion both in symplectic maps^ as well 
as in Hamiltonian flowsP^ 
, In order to compute the SALI for a given orbit one has to follow the time evolu- 

tion of the orbit itself and two deviation vectors which initially point in two different 



X 



directions. The evolution of these vectors is given by the variational equations for a 



flow and by the tangent map for a discrete-time system. At every time step the two 
vectors v\(t), V2(t) are normalized and the SALI is computed as: 



SALI(t) = min 



"It ^2 



Vllt v 2 



(1-1) 



where t is the continuous or the discrete time and II - II denotes the Euclidean norm. 
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The properties of time evolution of the SALI clearly distinguish between ordered 
and chaotic motion as follows: In the case of Hamiltonian flows or TV-dimensional 
symplectic maps with N ^ 2 the SALI fluctuates around a non-zero value for ordered 
orbits, while it tends to zero for chaotic orbits PS i n the case of 2D maps the SALI 
tends to zero both for ordered and chaotic orbits, following however completely 
different time rates, which again allows us to separate between these two cases alsoP 

We have recently begun to understand the different behaviors of the SALI in 
regions of order and chaos. In the latter case, we have been able to connect SALI's 
rapid convergence to zero, to the influence of the two largest positive Lyapunov 
exponents of the motion.^ In the present paper we shall study the behavior of the 
SALI in the case of ordered orbits. 



§2. The behavior of the SALI for ordered motion 



Let us try to understand why the SALI does not become zero in the case of 
ordered motion, by studying in detail the behavior of the deviation vectors. A 
suitable way to do this for conservative systems is to consider a non-trivial integrable 
Hamiltonian model whose orbits are bounded and lie on "nested" tori, which foliate 
all of the available phase space.™ 

An integrable such Hamiltonian system of 2 degrees of freedom possesses besides 
the Hamiltonian H a second independent integral F, in involution with H: 

{H,F} = 0, (2-1) 

where {•, •} denotes the usual Poisson bracket. In such systems, the motion lies in 
the intersection of both manifolds 

H = h, F = f, (2-2) 

where h, f are the constant values of the two integrals. Thus, the orbits in the 
4— dimensional phase space move instantaneously on a 2-dimensional "tangent" sub- 
space, which is 'perpendicular' to the vectors 

V-ff = (H x ,H y ,H Px ,H Py ), VF = (F x ,F y ,F Px ,F Py ), (2-3) 

x, y being the generalized coordinates of the system and p x , p y their conjugate 
momenta, while subscripts denote partial derivatives (e. g. H x = Qj^)- In fact, the 
motion may be thought of as governed by either one of the Hamiltonian vector fields 

Jh = (H Px ,H Py ,-H x , -H y ), fp = (F Px ,F Py ,-F x , —F y ). (2-4) 

The vectors VH, VF (and hence also fn, If) are linearly independent due to the 
functional independence of the two integrals at almost all points in phase space. So 
the corresponding unit vectors 

— Jh — If VH VF 

/ fl = -^-lV5, / F = 4^1VF, with VH = ^=^, VF=^^- (2-5) 

WM ||/f|| livtfll iivfii 
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can be used as a basis for the 4-dimensional space where the deviation vectors evolve. 
This basis is in general not orthogonal as 

(VH, VF) = (fn, Tf) = ^ + HyF lZ Hp d Px + HpyFpv (2-6) 

||ViT|| ||VF|| 

is not necessary zero. We note that ||Vi7|| = ||/h"||, ||VjP|| = and (•, •) denotes 
the usual inner product. Note also that from definitions (|2-3|) and (|2-4j) we get 

(VH, Th) = (VF, Tf) = 0, while (HI yields {VH, Tf) = (VF, Th) = 0. 

So, using vectors ()2-5j) as a basis for studying the evolution of a deviation vector 
Vi, we can write it as 

v\ = aifg + a 2 T + as^H + a 4 VF (2-7) 

with oi, a2, 03, «4 G K. The values of the coefficients ai, i = 1,2,3,4, at different 
times, give us a clear picture for the evolution of v\ . In the case of the 2D standard 
map for example, where ordered orbits lie on an invariant curve (ID torus), it has 
been shown both numerically and analytically 1 that any deviation vector (consid- 
ered as a linear combination of the vectors fjj, VJT using our notation), eventually 
becomes tangent to the invariant curve, tending to the tangential direction as n" 1 , 
with n being the number of iterations. 

Similarly, in the case of an integrable 2D Hamiltonian the deviation vector v\ 
tends to fall on the "tangent space" of the torus, spanned at each point by fn, /f, 
meaning that in Eq. 1)2- 7 Jl 03 — > 0, £14 — > 0, while the ai, 0,2 are, in general, different 
from zero. This is analogous to what has been found for the 2D standard map in 
Ref. |SJ) . As a model for studying this behavior let us consider the 2D Van der Waals 
Hamiltonian^ 

H(x,y,p x , Py ) = l -{p 2 x + P 2 y )-E{x 2 + y 2 ) + A{x % + y 6 ) + B(x 4 y 2 +x 2 y 4 ), (2-8) 

where E, A, B are real parameters. For B = 2>A and E G K the Hamiltonian (|2-8|) 
is completely integrable and the second integral of motion is given by® 

F(x,y,p x ,p y ) = (xp y -ypx) 2 . (2-9) 

In our calculations we consider the integrable case for A = 0.25, B = 2>A = 0.75 and 
E = -10" 8 . 

For different initial deviation vectors, we compute the time evolution of the 
coefficients ai, 0,2, 03, 04 of Eq. (|2-7jl . We find that in all cases 01, 02 remain different 
from zero, while 03, 04 tend to zero. A particular example is given in Fig. ^ By 
fitting the data of Fig. \T)p we see that 1 0.3 1 , |a-4 1 oc t . From Fig. ^ we conclude 
that any vector will eventually fall on the "tangent space" of the torus on which 
the orbit evolves. This "tangent space" is produced by vectors /#, fp, and so any 
deviation vector will eventually become a linear combination of these two vectors 
only. As there is no particular reason for two different initial deviation vectors to 
end up with the same values of a\, C12, the SALI (|1T|) will in general oscillate around 
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Fig. 1. The time evolution of the coefficients oi, £12, a3, 04, of the initial deviation vector with 
Oi = 1, d2 = 1, a-3 — 0, a4 = f (Eq. I2-7I 1. (a) oi (black line), ai (gray line), (b) |ei3 1 (black line), 
|<24j (gray line) in log-log scale. The orbit of the Hamiltonian 12-80 used, has initial condition 



£l 4 



= -0-6, y 
= 1.067- 



0, p x = 0, p y = 1.99416. 
-° " 5 that fit the data. 



In (b) we also plot the curves 



0.107 - t 



-0.<)«2 



a value different from zero. In other words the two vectors become tangent to the 
torus and fluctuate quasiperiodically about two different directions. This becomes 
evident in Fig. where we plot the time evolution of the SALI for an orbit with 
initial conditions x = —0.6, y = 0, p x = 0, p y = 1.99416 marked by a black point 
in the Poincare Surface of Section (PSS) of the system seen in Fig. [2Jd. The initial 
deviation vectors used are v{ = fn + If + VF (the time evolution of which is given 
in Fig. HJ) and vt = fn + If + Vff. 



§3. Conclusions 



In this paper we have analyzed the behavior of the SALI in regions of ordered 
motion, by studying the evolution of deviation vectors in the case of an integrable 
2D Hamiltonian system. Using a suitable basis of 4 vectors (|2-5|) we have shown 
that any pair of arbitrary deviation vectors tends to the tangential space of the 
torus, following a f 1 time evolution and having in general 2 different directions. 
This explains why for ordered orbits the SALI oscillates quasiperiodically about 
values that are different from zero. The same result is observed to hold in the 
case of ordered motion in a stability region of a non-integrable system;- where the 
presence of "islands" implies the existence of an additional approximate integral F, 
independent of the Hamiltonian. For Hamiltonian systems of more than 2 degrees 
of freedom we expect similar results. The only difference is that the "tangent space" 
is of higher dimension generated by the vectors fn , /fi > If 2 > • • with Fi, F2, . . . , 
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Fig. 2. The time evolution of the SALI (a) for the ordered orbit marked by a black point in the 
PSS y = of the system (b). 

being the additional (approximate or not) integrals of the motion. 
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